*Manuscript 

Click here to view linked References 


The Structure of Titan’s Atmosphere from 
Cassini Radio Occultations 

Paul J. Schinder a , 

^Center for Radiophysics and Space Research, Cornell University, Ithaca, NY 

14853 (U.S.A.) 

F. Michael Flasar b 

h NASA Goddard Space Flight Center, Code 693, Greenbelt, MD 20770 (U.S.A.) 

Essam A. Marouf c 

c Department of Electrical Engineering, San Jose State University, San Jose, CA 

95192 (U.S.A.) 

Richard G. French d 

d Department of Astronomy, Wellesley College, Wellesley, MA 02481 (U.S.A.) 

Colleen A. McGhee 6 

e Department of Astronomy, Wellesley College, Wellesley, MA 02481 (U.S.A.) 

Arvydas J. Kliore f 

f Jet Propulsion Laboratory, Pasadena, CA 

Nicole J. Rappaport 8 

Preprint submitted to Icarus July 12, 2011 


s Jet Propulsion Laboratory, Pasadena, CA 


Elias Barbinis 11 

11 Jet Propulsion Laboratory, Pasadena, CA 


Don Fleischman 1 

1 Jet Propulsion Laboratory, Pasadena, CA 


Aseel AnabtawB 

J Jet Propulsion Laboratory, Pasadena, CA 


Number of pages: 38 
Number of tables: 3 
Number of figures: 12 

Email addresses: schinder@astro.cornell.edu (Paul J. Schinder), 
f.m.flasar@nasa.gov (F. Michael Flasar), emarouf@email.sjsu.edu (Essam A. 
Marouf), rfrench@wellesley.edu (Richard G. French), 
cmcghee@firstclass.wellesley.edu (Colleen A. McGhee), 

akliore@jpl.nasa.gov (Arvydas J. Kliore), Nicole.J.Rappaport@jpl.nasa.gov 
(Nicole J. Rappaport), elias.barbinis@jpl.nasa.gov (Elias Barbinis), 
don.u.fleischman@jpl.nasa.gov (Don Fleischman), 
Aseel.Anabtawi@jpl.nasa.gov (Aseel Anabtawi). 


2 




Proposed Running Head: 

The Structure of Titan’s Atmosphere 


Please send Editorial Correspondence to: 


Paul J. Schinder 
1034 Karen St. 
Boalsburg, PA 16827 


Email: schinder@astro.cornell.edu 
Phone: (814) 466-2150 


3 



ABSTRACT 


We present results from the two radio occupations of the Cassini spacecraft 
by Titan in 2006, which probed mid-southern latitudes. Three of the ingress 
and egress soundings occurred within a narrow latitude range, 31-34 ° S near 
the surface, and the fourth at 52.8° S. Temperature - altitude profiles for all 
four occupation soundings are presented, and compared with the results of the 
Voyager 1 radio occupation (Lindal et ah, 1983), the HASI instrument on the 
Huygens descent probe (Fulchignoni et ah, 2005), and Cassini CIRS results 
(Flasar et ah, 2005; Achterberg et al., 2008b). Sources of error in the retrieved 
temperature - altitude profiles are also discussed, and a major contribution is 
from spacecraft velocity errors in the reconstructed ephemeris. These can be 
reduced by using CIRS data at 300 km to make along-track adjustments of the 
spacecraft timing. The occupation soundings indicate that the temperatures 
just above the surface at 31-34° S are about 93 K, while that at 53° S is about 
1 K colder. At the tropopause, the temperatures at the lower latitudes are all 
about 70 K, while the 53° S profile is again 1 K colder. The temperature lapse 
rate in the lowest 2 km for the two ingress (dawn) profiles at 31 and 33° S 
lie along a dry adiabat except within ~ 200m of the surface, where a small 
stable inversion occurs. This could be explained by turbulent mixing with low 
viscosity near the surface. The egress profile near 34° S shows a more complex 
structure in the lowest 2 km, while the egress profile at 53° S is more stable. 


Keywords: TITAN; ATMOSPHERES, STRUCTURE 
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1 Introduction 


Titan is unique among the moons of the solar system in having a thick, ex- 
tended atmosphere, with a surface pressure exceeding 1 bar (1000 hPa). N 2 
is the dominant constituent, followed by CH 4 , which has a mole fraction of 
0.015 in the stratosphere increasing at lower altitudes in the troposphere to 
several percent at the surface (Niemann et al., 2005, 2010). The most de- 
tailed measurements of temperatures in Titan’s troposphere and middle at- 
mosphere (stratosphere and mesosphere) have been made by instruments on 
the Cassini orbiter and Huygens probe. The Huygens Atmosphere Structure 
Instrument (HASI) obtained temperatures along the descent trajectory near 
10° S (Fulchignoni et ah, 2005). Between 1400 km and 160 km, atmospheric 
mass density was calculated from accelerometer data and using a drag co- 
efficient characterizing the probe. Application of the barometric law yielded 
pressure, and the ideal gas law provided temperature. Below 160 km, with 
the parachutes deployed, temperature and pressure were directly measured by 
sensors. Vertical resolution was determined by the time sampling, and it was 
20 km between 400 km and 1400 km altitude, decreasing to 1 km near 160 
km (Fulchignoni et ah, 2005). At lower levels it was finer, approximately 10 
m just above the surface (Tokano et ah, 2006). 

In the thermal infrared, the Composite Infrared Spectrometer (CIRS) on the 
Cassini orbiter has mapped temperatures globally (Flasar et ah, 2005; Achter- 
berg et ah, 2008b). In principle, it can retrieve temperatures between 300 mbar 
(~ 28km) in the upper troposphere (Flasar et ah, 2004) and ~ 3/ibar (~ 470 
km, Achterberg et ah, 2008b). To date, however, little has been reported on 
temperatures in the lower stratosphere between 100 and 10 mbar (50-100 km). 
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This is because, in addition to collision- induced absorption from N 2 — N 2 and 
N 2 — CH4 pairs, aerosols and condensates also contribute to the far-infrared 
opacity, and the latter have not been fully analyzed yet. Because of the low- 
pass filtering in the radiative transfer process and the distance of the spacecraft 
from Titan when limb-viewing observations are made, the vertical resolution 
in the CIRS temperature retrievals is typically a pressure scale height, e.g., 
~ 40 km in the middle atmosphere. Temperatures at higher altitudes in the 
mesosphere (~ 500 km) have also been retrieved from solar and stellar oc- 
cultations by the ultraviolet spectrometer on Cassini (e.g., Shemansky et ah, 
2005) and above 400 km from stellar occultations observed from Earth-based 
observatories (e.g., Hubbard et al., 1993; Sicardy et al., 2006). 

Compared to the other Cassini orbiter experiments, radio occultations are 
the only means of obtaining temperatures in the lowest 1.5 scale heights of 
Titan’s atmosphere, as well as the lower stratosphere at present. Moreover, 
the retrieved temperatures have much higher vertical resolution, which can, 
for example, resolve thin layers and waves with small vertical wavelengths. 
The technique, in which a radio signal from a spacecraft is refracted by the 
atmosphere of a target body before being recorded at a station on the Earth, 
can determine the structure of an atmosphere over several decades of pressure. 
It has been successfully used for many bodies in the solar system, including 
Jupiter and its satellites (Eshleman, 1975; Hinson et al., 1997), Saturn (Lindal 
et al., 1985), Uranus (Lindal et al., 1987), Neptune and Triton (Lindal, 1992), 
Mars (Cahoy et al., 2006), and Venus (Jenkins et al., 1994; Hinson and Jenkins, 
1995). Prior to Cassini, the Voyager 1 spacecraft was occulted by Titan in 
1979, early in Titan’s northern spring (Lindal et al., 1983). The two raclio- 
occultation soundings were both near the equator (Table 1). 
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[Table 1 about here] 


In 2006 (during Titan’s northern winter), the second and third radio occul- 
tations of a spacecraft by Titan occurred. The Cassini spacecraft, in orbit 
around Saturn, was occulted twice by Titan, and Cassini obtained soundings 
by radio occultation at more southerly latitudes (Table 1) than Voyager. The 
first occultation occurred on orbit 22 during the T12 flyby on March 19. (Tn 
designates the sequence of the targeted Titan encounters. In order these are 
TA, TB, TC [when the Huygens probe descended into Titan’s atmosphere], 
T3, T4, T5, etc.) The second occultation, T14 on orbit 24, occurred on May 
20. The Cassini radio-occultation experiment differs in two important ways 
from that on Voyager (Kliore et ah, 2004). First, for the first time, there are 
three frequencies: S-band (2.3 GHz), X-band (8.4 GHz), and Ka-band (32 
GHz); Voyager had only S and X. Second, the Cassini Ultra Stable Oscillator 
(USO) is a factor ~ 25 more stable than Voyager’s. The latter permits sound- 
ing the neutral atmosphere higher than the 200 km achieved with the Voyager 
occultation soundings (Lindal et ah, 1983). Since error in the bending angle 
(the angle between the initial ray direction as it enters the atmosphere and 
the final direction as it leaves) is directly related to error in the frequency, and 
the bending angle is smaller at high altitudes, a more stable frequency source 
with smaller errors gives us the capability to retrieve the temperature profile 
at higher altitudes. 

Here we report the results of the analysis of the 2006 occupations. Two addi- 
tional occupations at higher latitudes in both hemispheres occurred in 2007 
during Cassini’s primary mission. During the extended mission, one occupa- 
tion at mid latitudes occurred in November 2008, and another ingress only 
occupation at high northern latitude in June, 2009. These will be the subject 
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of a separate paper. 


2 The 2006 occultations 

2.1 Procedure and Results 

During the T12 occultation, Cassini transmitted at three frequencies refer- 
enced to the on board USO (Kliore et al., 2004). The signal was collected 
using three antennas at the Deep Space Network station at Goldstone, Cali- 
fornia (DSS-14, 25, and 26), and two antennas at the DSN station at Madrid, 
Spain (DSS-55 and 63). Data were simultaneously recorded at the S (2.3 GHz), 
X (8.4 GHz), and Ka (32 GHz) bands. The T14 occultation was also recorded 
at the three bands by the Madrid Deep Space Network station over the 70-m 
antenna DSS-63 (S and X) and the 34-m DSS-55 (X and Ka). During both oc- 
cultations the boresight of the high-gain antenna on Cassini was continuously 
pointed at the predicted position of the virtual image of the Earth, that point 
on Titan where the refraction of the atmosphere will bend the ray so that it 
will intercept the Earth. At reception, the received signal from the spacecraft 
was “steered” by mixing with a reference frequency, which was computed by 
using a model Titan atmosphere and a code that predicts the received fre- 
quency of rays passing through that atmosphere, into a low frequency band 
that can be digitally recorded. The resultant mixed signal was then sampled 
at three different rates, 1, 16, and 50 kHz. 

To determine the structure of the atmosphere of Titan with these data, we 
begin by Fourier analyzing the samples within a time window to determine the 
frequency at that time. Throughout this paper, we use the data sampled at 



16 kHz, and use 4096 data points (.256 seconds) in a sliding window, with the 
window sliding 512 points (.032 seconds) between each frequency determina- 
tion. Adding this to the reference frequency allows us to reconstruct the time 
series of the sky frequency, the original frequency of the signal as received at 
the DSN antenna. By averaging a series of those rays that don’t enter Titan’s 
atmosphere (“the baseline”), prior to the beginning of the actual sounding by 
the atmosphere (or after it, in the case of an egress sounding), and knowing the 
positions and relative velocity velocities of the spacecraft and DSN receiving 
station, we can determine the frequency at which the spacecraft transmitted 
in its rest frame. The fractional difference between the transmitted and the 
sky frequencies provides the initial ray direction, or, more precisely, the an- 
gle between the ray direction and the spacecraft velocity vector. Knowledge 
of the initial ray direction and final ray direction allows us to determine the 
refractivity profile of the atmosphere, as explained below. 

We assume that the atmosphere of Titan is windless and spherically sym- 
metric. Under these conditions, the Abel transform technique (Phinney and 
Anderson, 1968; Fjcldbo et al., 1971) can be used to retrieve the refractiv- 
ity structure of the atmosphere. (The refractivity N is related to the index 
of refraction n by N = 10 6 (n — 1).) We can alternatively use ray tracing, 
a more general technique that assumes the atmosphere is axisymmetric and 
barotropic. Ray tracing provides additional useful information, such as the 
actual paths of the rays through the atmosphere, but at a much higher com- 
putational cost. However, the refractivity structure found by our ray tracing 
and Abel codes are identical for a spherically symmetric atmosphere, and the 
results presented below were obtained with our Abel code. 

Figure 1 depicts the fractional power transmitted through Titan’s atmosphere 


9 



during the ingress and egress of T12 and T14. Atmospheric absorption, in 
addition to refractive defocusing, occurs at the X- and Ka bands. The signal 
was tracked all the way to the surface for S and X bands, while the Ka signal 
disappeared at about 10 km above the surface. The primary contribution to 
absorption is from pressure-induced absorption by N 2 (Ho et ah, 1968; Bo- 
rysow and Frommhold, 1986; Lindal et ah, 1985). A more detailed analysis of 
atmospheric absorption on Titan at radio wavelengths will be given elsewhere. 

[Figure 1 about here.] 

The Abel (or the ray-tracing) code derives, from the top down, a profile of 
the refraetivity of the gas as a function of altitude. We obtain the number 
density by dividing by the refraetivity per molecule, and the mass density by 
multiplying by the mean molecular weight of the gas. We assume the atmo- 
sphere of Titan is mostly N 2 , with a small admixture of CH 4 at a level of 
1.4 x 1CH 2 by number (Niemann et ah, 2005). The corresponding refraetivity 
per molecule used for N 2 is 1.1 x 10 -17 cm 3 per molecule, and that for CH 4 is 
1.63 x 10~ 17 cm 3 (Achtermann et ah, 1991; Orcutt and Cole, 1967; Essen and 
Froorne, 1951). The equation of hydrostatic equilibrium is integrated down- 
wards from a chosen starting altitude and temperature. We use the Bender 
equation of state, which takes into account the departure of the equation of 
state from a perfect gas at the low temperatures and high pressures character- 
istic of the lower troposphere and conveniently permits the treatment of gas 
mixtures (Bender, 1971, 1973; Buhner et ah, 1981). 

[Figure 2 about here.] 

Figure 2 shows the latitude and longitude of the ray periapses as they descend 
toward (or ascend from, during egress) the surface of Titan. The soundings 
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usually cover several degrees of latitude and longitude as the virtual image 
of the spacecraft as seen from the earth descends toward (ascends from) the 
surface. Figure 3 shows the temperature-altitude profile for the region between 
0 and 300 km for the T12 and T14 soundings. The integration was started as 
near to 300 km as possible at 185 K; at that altitude the CIRS profiles we 
use to obtain a starting value have an error of about ±1K (Achterberg et ah, 
2008b). Appendix A contains tables that give the values we obtain for the 
pressure, temperature, and refraetivity at selected altitude levels for all four 
soundings. As described below, we have adjusted the spacecraft ephemeris so 
that the temperature profile in the stratosphere closely matches that obtained 
by CIRS observations. 


[Figure 3 about here.] 

The profiles given in the tables and the figures have been retrieved from 16 kHz 
X-band data. We have found that profiles in the neutral atmosphere below 300 
km retrieved from X and Ka-band data are little affected by refraction from 
electrons in Titan’s weak ionosphere. The S-band data, however, are noticeably 
affected by the presence of electrons. It is well established that this effect can be 
removed by using two frequencies. For example, the refraetivity derived from S 
band data ris is assumed to consist of two parts, no independent of frequency 
and rii(f) a function of frequency, and similarly for the X band refraetivity rix- 
It follows that n 0 = (n s - M/s)M(/x)]nx)/( 1 - M/s)M(/x)]) (Schunk 
and Nagy, 2000). Since n t oc /~ 2 for free electrons, we can use the X and S 
band, or any other two bands, to extract the refraetivity no of the neutral 
gas alone. We have found that the retrieved temperatures using the X- or 
Ka-band (above 10 km) alone are virtually identical to the profiles obtained 
by combining X and S bands in the altitude range we present here. 
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2.1.1 Diffraction 


Because we use geometric optics in the retrievals, diffraction limits the vertical 
resolution of structure in the atmosphere to the diameter of the first Fresnel 
zone (Fjeldbo and Eshleman, 1969). In the stratosphere, the diameter is close 
to the free-spaee value, 2 V\D, where A is the wavelength of the radio signal 
and D the distance of Cassini from Titan’s occulting limb. For X band, this is 
typically ~1 km for the T12 and T14 occupations. Differential refraction in a 
neutral atmosphere reduces the vertical scale by a factor rs_/ 1/V(1 + Da/H), 
where H is the pressure scale height and a is the bending angle between 
the incoming and outgoing ray asymptotes (Haugstad, 1978). In the lower 
troposphere Da/H » 1, and the vertical diameter of the first Fresnel zone 
is ~ 2 \J XH / a r \ j 300 m for X band, independent of D. 

Titan has a solid surface, and the surface can cause diffraction in the last 
~ 300 m of the temperature profile. Since diffraction would cause a devia- 
tion from the geometric optics solution, the Abel transform (or ray tracing), 
which assumes geometric optics, might misinterpret diffraction as a temper- 
ature structure right at the surface. Since we use exclusively the frequency 
record and not the power record, we are sensitive only to the deviation in 
the actual recorded frequency from the geometric optics prediction. It should 
be noted that there is no convincing signal in either the power record or the 
frequency record that diffraction is occurring at the surface, but this may be 
due to the fact that the power is also dropping rapidly at the surface due to 
refractive defocusing (for the X and Ka bands, pressure induced absorption by 
N 2 is also a factor). Anything, not just diffraction, that could cause a devia- 
tion from the geometric optics solution in a spherically symmetric atmosphere 
(such as we are assuming for the Abel transform) could cause deviations of 
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the computed temperature profile from the actual temperature profile. These 
include actual atmospheric structure that is not spherically symmetric, both 
large scale zonal and longitudinal deviations from spherical symmetry and 
small scale deviations (waves, etc.). The reverse is also true: structure that 
actually exists in the atmosphere near the surface (local surface inversions, for 
example, or deviations caused by terrain) may be misinterpreted as the effects 
of diffraction. 

French and Gierasch (1976) show that the effects of diffraction are greatly 
reduced in the presence of a thick atmosphere, such as found on Titan. One 
way to remove the effects of diffraction from the recorded data and therefore 
be able to distinguish between real surface effects and diffraction is backprop- 
agation (Karayel and Hinson, 1997; Marouf et al., 1986). Backpropagation in 
effect synthetically shifts the trajectory of the spacecraft to one closer to the 
target body so that the Fresnel zone becomes arbitrarily small. In this case, 
backpropagation must be done taking into account the presence of a thick 
atmosphere between the DSN station recording the frequency data and the 
spacecraft. We have devised a procedure to do this, but the details and the 
results of both backpropagation and the results of forward calculations using 
a simplified geometry will be held for a future paper to allow a full discussion 
of the problem. The results we present here are not backpropagated and the 
vertical resolution is diffraction limited. 


2.2 Sources of error 

In this section, we will discuss some of the sources of error that affect the 
profiles, particularly those that affect the profiles preferentially above 70 km, 
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and how we have corrected for them when possible. 


2.2.1 Hydrostatic equilibrium starting point 

First, the computation of the temperature and pressure profiles from the re- 
fractivity profiles starts at a given altitude with an arbitrarily chosen temper- 
ature. 


[Figure 4 about here.] 

Figure 4 shows the effect of starting the hydrostatic equilibrium calculation 
at 300 km at three different temperatures, 175, 185, and 195 K, for T14 egress 
using X band data. Note that the shapes of the profiles remain the same, 
and the divergence decreases (approximately as 1/P, see Lipa and Tyler, 
1979) as we descend, becoming negligible at about 150 km. Figure 4 also 
shows the effect of starting the hydrostatic equilibrium calculation at 148 
km at 172.5 K, approximately the same as the first direct measurement by 
the HAS1 instrument on Huygens, again for T14 egress. Both profiles agree 
well, although the radio occultation profile does not show the abundant small 
scale structure seen in the Huygens profile. This is not surprising, because the 
Huygens in situ measurements detect local structure, whereas the occultation 
measurements detect an average along the ray path in a shell of thickness 
comparable to the vertical scale of the structure being resolved (provided it is 
larger than the Fresnel scale). This is roughly the geometric mean of Titan’s 
radius (Rt = 2575 km) and the vertical scale. In Figure 4 the occultation 
retrievals would effectively average structure observed by Huygens on 10-km 
scales over a distance ~ \/10 x 2575 ~ 160 km. If the horizontal scale of the 
small-scale structure measured by Huygens is much smaller, then it will be 
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smoothed out. 


In order to minimize the arbitrary starting point as a possible source of error 
in our profiles, we use CIRS profiles (Achterberg et ah, 2008b) to set the 
temperature at 185 K at an altitude of 300 km. At that altitude the CIRS 
retrievals have an error (from instrument noise) of approximately ±1K. 

2.2.2 Baseline 

The start of every ingress occultation and end of every egress occultation 
includes a baseline , a period of time in which the signal from the spacecraft is 
tracked in “free space”. This is done so that the spacecraft USO frequency in 
its rest frame can be determined from the data, to enable subsequent analysis. 
However, there are choices to be made in how to handle the baseline, and 
what to do with the inevitable deviations from the ideal of transmission in a 
vacuum. 

The choice begins with an examination of the Doppler residual frequency, the 
difference between the actual earth received frequency and the earth received 
frequency predicted if the path between the spacecraft and the DSN antenna 
were in vacuum. For this purpose, the spacecraft oscillator frequency is deter- 
mined from the first point in the time series of earth received frequencies vs. 
time. 


[Figure 5 about here.] 

Figure 5 shows an example for the T14 ingress occultation. There are several 
features to note. First, the Doppler residuals don’t lie exactly on zero, as they 
would if both the receiver and transmitter are in vacuum. Instead, they show 
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1) a subtle linear trend away from zero, 2) subtle small scale structure, and 3) 
a narrow spread in frequency due to thermal noise. Of course, the path from 
spacecraft to receiver isn’t entirely in vacuum (the signal passes through the 
Earth’s ionosphere and atmosphere, for example), so deviations from the ideal 
should be expected. 

We next choose an interval in time. We usually try to choose an interval as long 
as possible that doesn’t have obvious problems (excess curvature, discontinu- 
ities, etc.). We also need to avoid the target’s ionosphere, if any. A baseline of 
more than 200 s is preferred, although not always possible. In the case shown 
in Fig. 5, we choose the interval 48000-48400 s. For some soundings, such as 
T12 ingress, the baseline is short for reasons beyond our control, and we have 
to accept the greater (mostly stratospheric) uncertainty that results. The er- 
rors caused by any sort of frequency uncertainty are always larger where the 
density and the bending angle are smaller, so they affect high altitudes the 
most. 

Next, we numerically integrate the Doppler residuals to derive a phase, and 
then fit the phase with a cubic polynomial using a linear least squares fit. The 
derivative of this cubic gives us a quadratic fit to the residual frequency. We 
then determine the spacecraft oscillator frequency from this fit. Ordinarily we 
then remove the linear drift using the linear part of the fit, although occasion- 
ally the quadratic part is big enough so that it must be removed as well. In 
this way we get as close as possible to the ideal of the Doppler residuals being 
zero outside the atmosphere. 

However, the baseline and what we do with the baseline is a subjective choice 
to some degree, and that choice and how we correct the data will affect the 
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subsequent temperature - altitude profile. In order to minimize the effects of 
both the baseline choice and the spacecraft ephemeris errors which we describe 
next (which also cause errors in frequency), we choose to use CIRS profiles 
(Achterberg et ah, 2008b) and adjust the ephemeris by introducing a timing 
offset between the time a given ray leave the spacecraft and the time at which 
we query the ephemeris for the spacecraft’s position and velocity for that 
ray. This allows us to bring our stratospheric temperature - pressure profile 
into agreement with the CIRS profiles, while eliminating both the systematic 
errors caused by the choice of the hydrostatic equilibrium starting point and 
the choice of the baseline, leaving only random errors. 


2.2.3 Ephemeris errors 

The ephemerides provided by JPL come in the form of binary kernels which 
contain fits to the actual trajectory of the spacecraft, Saturn, and Titan for 
a certain time interval. The errors in both the position and the velocity of 
the spacecraft change with time along the trajectory. The position errors for 
the “reconstructed” kernels (kernels which contain the past trajectory of the 
spacecraft based on the best available tracking data) are on the order tens of 
meters near Titan. The velocity errors are on the order of a few mm/s. 

We have equipped our code with the ability to do an along-track “time shift” ; 
in other words, we query the ephemeris at time t + At and use that state 
vector as if it were the correct one for time t. We perform this shift only for 
the spacecraft state vector, and not for Saturn, Titan, or the DSN antenna. 
This is a simple, one parameter attempt to estimate the effect of ephemeris 
errors. More elaborate schemes can be envisioned (including, for example, 
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using information from the atmospheric occultation to improve the ephemeris), 
but are much harder to implement. 

[Figure 6 about here.] 

Figure 6 shows the temperature profiles obtained for various values of At. The 
effects of even a few millisecond time shifts are quite apparent above 70 km, 
with large temperature differences above 150 km. The rejoining of the profiles 
at 300 km at 185 K is an artifact of the fact that the hydrostatic equilibrium 
calculation was started there. 

Using a time shift means that both the position and the velocity of the space- 
craft change from the nominal value. It turns out that the velocity change 
is the controlling factor. We have tested this by suppressing the change in 
either the position or the velocity caused by the time shift. While the shift in 
position of a few tens of meters causes only a very small (< 1 K) temperature 
shift at high altitudes, the shift in velocity of ~ 1 mm s” 1 causes almost all of 
the shift shown in Fig. 6. 


[Figure 7 about here.] 

We have obtained from the Cassini Navigation (NAV) team a set of covariance 
matrices describing the error in position and velocity of the reconstructed SP 
(Spacecraft-Planetary) kernel from JPL’s Navigation and Ancillary Informa- 
tion Facility (NAIF) for both of these occupations. With this, we are able 
to make a more realistic appraisal of the effect of ephemeris errors on the 
stratospheric temperature profile. 

Shifting the velocity affects the temperature profile in two ways. First, the 
spacecraft USO frequency is determined from the sky frequencies recorded on 
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the ground. The relation between the spacecraft oscillator frequency f s and 
the sky frequency f a in the rest frame of the DSN antenna is 


fa _ 1 

fs 1 — (v s • n)/c 

to first order in v/c ., where n is the unit vector along the initial direction of 
the ray starting at the spacecraft and ending at the ground antenna, v s is the 
spacecraft velocity, and c the speed of light. As described above, we essentially 
average this equation over the baseline to get the constant spacecraft frequency 
f s . So if the ephemeris returns a velocity of v + 5v when the spacecraft velocity 
is actually v, the spacecraft oscillator frequency will be slightly in error by an 
amount 


fs = fs 0 -<fa(Svn)/c> (2) 

where f so is the actual spacecraft frequency and the angle brackets denote the 
averaging procedure. 

Once the spacecraft oscillator frequency is obtained, we use it and the recorded 
ground frequencies to find the angle between the initial ray direction and the 
spacecraft velocity vector. The angle can easily be found from 

"A fa fs /g\ 

C fa 

The ray then must be bent by just the right amount to arrive at the DSN 
antenna with the measured frequency. This means that if the ray direction is 
slightly off from its actual direction due to ephemeris errors, the structure of 
the atmosphere computed by using the ray must change slightly to bend it in 
the right direction. This has a direct effect on the refractivity profile and an 
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indirect effect on the temperature profile. 


To begin, we choose a fixed direction at random. We then multiply that direc- 
tion by a changing magnitude calculated from the covariance matrix. Because 
the covariance matrix changes with time, changing the magnitude of the veloc- 
ity error gives a more realistic appraisal of the error caused by the ephemeris. 
The fixed direction for the velocity error (but not for the velocity itself) isn’t 
as realistic, but the actual change in direction of the total velocity is con- 
trolled by celestial mechanics, and a fixed direction for the error vector is a 
more reasonable approximation (for error analysis only) for the short dura- 
tion of an atmospheric occultation than to allow, say, the direction of the error 
component to randomly change with time. 

The covariance matrices are 1 er errors, so we performed a Monte Carlo esti- 
mation of the errors in the temperature - altitude profile by performing 100 
cases each of 3 a velocity errors. We again started the hydrostatic equilibrium 
integration at 300 km at 185 K. Figure 7 shows the resulting temperature 
profiles for T14 egress in X band using 16 kHz data. With 3 a, the maximum 
error is around 20 K. Once again, the joining of the profiles at 300 km is an 
artifact of the initial conditions. 

The la velocity errors provided by the NAV team for the T14 reconstructed SP 
kernel were 2-4 mm s^ 1 during the course of the occultation, while for T12 
the errors were 0.5 - 1 mm s -1 . Note that an error ~ 1 mm s~ x corresponds 
to an additional fractional frequency shift of 0.1 cm s -1 / 3 x 10 10 cm s -1 ~ 
3 x 10 -12 , or roughly 25 mHz for X band. In the lower stratosphere and tropo- 
sphere the fractional shift in frequency from refraction is much larger, ~ 10 - ', 
so the effect of the velocity error damps out in this region, as Fig. 7 shows. 
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Timing errors like the ones we discussed above are possible if the clock that 
timestamps the recordings in the radio science receivers are not in sync with 
the ephemeris time in the SP kernel. This may be possible, but the computers 
providing the timestamps for the radio science data hies are kept in sync with 
a local highly accurate time reference. Our most likely source of this class of 
error for determining the stratospheric temperature profile is timing errors in 
the ephemeris itself. 


[Figure 8 about here.] 

We have attempted to correct for the errors introduced by the velocity shifts 
by making use of retrieved CIRS temperatures near 300 km (Achterberg et al., 
2008b). To make the approach tractable, we do not alter the spacecraft veloc- 
ities in the ephemeris, but instead introduce a single time shift in the ingress 
and egress portions of each occupation. We adjust the time shifts until the 
difference between our profiles and the CIRS profiles are small, bringing our 
stratospheric temperature - pressure profiles into close agreement with the 
profiles from CIRS. Figure 8 illustrates the quality of the fitted profiles. It 
is possible with along-traek shifts within the quoted errors of the spacecraft 
ephemeris to align the RSS temperature profiles with the CIRS profiles to 
good accuracy. 

The CIRS temperature profiles used are zonal averages based on a combina- 
tion of limb and nadir observations from July 2004 to May 2007 (Achterberg 
et ah, 2008a, b). The limb profiles are limited and sample longitude poorly, 
but pressures 0.2 to 5 mbar are dominated by nadir-sounding observations, 
which sampled longitude very well and typically consist of averages of ~ 100 
or more retrieved temperature profiles within 5°-latitude bins. This pressure 
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range dominates that of the fitting procedure illustrated in Fig. 8, and the er- 
rors in temperature over these pressures from the CIRS instrument noise are 
approximately 0.5 K for a single retrieved profile. When zonally averaged, this 
reduces by a factor of 10 or more. The radio-occultation retrievals are not zonal 
averages, but sounded over a narrow range of longitudes (Fig. 2). The CIRS 
temperatures near 1 mbar, however, indicate small variations around latitude 
circles at the latitudes of the occupation soundings, with a standard devia- 
tion of < 0.6K. Temporal variations in the CIRS profiles over the three-year 
period also appear to be small, with seasonal variations becoming apparent 
in 2008-2009. This is understandable, because the earlier observations shortly 
after the northern winter solstice, when there is an extremum in the solar 
declination. 


2.2.4 Random noise in the frequency 

Figure 5 shows a spread of frequencies around Doppler zero, in this case with 
a width of ~ 40 mHz. The cause of this random error is thermal noise in the 
system. 

To a certain extent the size of this random component is under our control. 
Fourier transforming more data for each frequency point improves the signal 
to noise ratio, but including too much data creates artifacts that can affect 
the profiles. 

To estimate the effect that the random noise has on our profiles, we add white 
noise to the frequency record with the same ~ 40 mHz width, effectively 
doubling the noise, and do 100 trials. 
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[Figure 9 about here. 


The results are shown in Fig. 9. For T14, the errors caused by white noise are 
roughly half as large as the errors caused by ephemeris errors (see also Fig. 
3). 

2.2.5 Other sources of error 

Other sources of error include the assumed composition of the atmosphere and 
the value of the refractivity per molecule that we use. If some other constituent 
of the Titan atmosphere that we didn’t consider had an unusually large re- 
fractivity per molecule, then we would overestimate the density (the number 
of molecules required to obtain the observed refractivity), and this would af- 
fect our temperature profiles. However, the temperature profiles determined 
by radio occultation and other techniques agree well in the lower atmosphere, 
so errors of this sort are probably not large. As we showed above, the starting 
value of the hydrostatic equilibrium calculation can cause substantial errors 
in the stratosphere. However, the structure of the atmosphere at 300 km is 
well determined by other techniques, and we are using one of them to set the 
temperature at 300 km, so in practice this should not be a significant source 
of error. 


3 Discussion 

The left panels (a, c, e) of Fig. 10 display the T12 and T14 profiles together. 
The most notable feature is how similar they are, although differences do exist. 
Above 100 km, T14 egress is roughly 2.5 K warmer than the profiles at T14 
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ingress and T12 ingress, even though all are at nearly the same latitude. This 
is within the l-er errors of these retrievals from thermal noise at these altitudes 
(Fig. 3). Nadir sounding by CIRS (Achterberg et al., 2008a) indicates little 
zonal variation near 1 mbar at these latitudes. However, it is possible that the 
fits to the CIRS profiles in the upper stratosphere by simply shifting the tim- 
ing of the spacecraft ephemeris has introduced some residual systematic error. 
The right panels (b, d, f) compare the T14 retrievals with those from Voyager 
radio occultations and HASI. The Voyager profiles, although equatorial and 
at a different season, are nonetheless close to those from T12 and T14. The 
HASI profile is consistently warmer in the stratosphere, although this may in 
part be attributable to adjusting the spacecraft ephemeris and boundary con- 
dition to the CIRS profiles in inverting the radio-occultation data. There is a 
discrepancy between HASI and the CIRS retrievals at low latitudes that is not 
understood. The HASI profiles are ~ 10 K warmer in the upper troposphere, 
and calculated synthetic spectra using these profiles produce infrared spectra 
that have radiances near 7 pm that are much too high compared to CIRS 
observations (Coustenis et ah, 2005; Ferri et ah, 2006; Flasar and Achterberg, 
2009). 

The three profiles at 32-34° S are virtually identical in the tropopause region 
(panel c, 30-55 km). The profile at the higher latitude, 53° S, is approximately 
0.8 K cooler near 40 km. The 32-34° S profiles very close to those from HASI 
at these altitudes, but up to 0.4-0. 8 K cooler than those from Voyager radio 
occultations (panel d). 

Closer to the surface (panel e), the profiles at 32-34° S are close to dry adi- 
abatic in the lowest 2 km, and gradually becoming more stable above that. 
This can be seen in Fig. 11, which depicts the vertical gradients of temperature 
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in the lowest 10 km. In the T12 ingress and T14 ingress profiles, there is an 
indication of more stable behavior below the adiabatic layers just above the 
surface. Both ingress soundings occurred near the dawn terminator, and this 
stability may reflect surface cooling overnight. A stable nocturnal inversion 
would not be radiatively driven, but more likely induced by turbulent mixing 
with low effective viscosity (see Appendix B). An adiabatic layer extending 
up to 2 km is more or less consistent with the Voyager results of Lindal et al. 
(1983), but it is much greater than the 300 m height inferred from the HASI 
data near 10° S (Tokano et ah, 2006). However, one should bear in mind that 
the HASI measurements are local, whereas the radio occultation temperature 
at a given altitude is really an average of the atmospheric structure over the 
ray path, as discussed earlier in section 2.2.1. 

The 32-34° S temperatures and those from HASI and the Voyager occultation 
are virtually identical at the surface (panel /). At higher altitudes, the HASI 
profile is relatively warmer, about 0.7 K at 6 km. The Voyager profiles are not 
as warm, and lie closer to the Cassini radio-occultation temperatures. The 
profile at 53° S, near the evening terminator, is more stably stratified below 
3 km than the lower- latitude Cassini profiles. At 0 km, relative to a 2575-km 
sphere, it is 0.9 K cooler than the other three (panel e). 

Over much of the region depicted in Fig. 11, the temperature gradients ex- 
hibit a fine-scale “spiky” structure. This fine structure may reflect ray crossing 
attributable to small-scale vertical variations in the refractivity. French and 
Lovelace (1983) have examined the distortion that a sinusoidal variation in 
refractivity of sufficiently large amplitude can produce in temperature profiles 
retrieved from stellar occultations using geometric optics. Appendix C briefly 
describes the effect for for a sinusoidal variation added to an isothermal at- 


25 



mosphere. A real atmosphere will probably have a spectrum of such waves, 
and a more complex structure. We have made a preliminary study using the 
backpropagation procedure mentioned in section 2 . 1 . 1 , and we find that the 
spikes have smaller amplitudes as the observer’s trajectory is synthetically 
shifted toward the occulting limb, which suggests that ray crossing is indeed 
occurring. Note that the spiky behavior in Fig. 11 is largely absent from the 
adiabatic layers. This may be because the more rapid vertical mixing in these 
layers has smoothed out the small perturbations. A sharp feature at 1 km al- 
titude appears to be “embedded” within the adiabatic layer of the T14 egress 
profile, which is near the evening terminator. It might represent a thin stable 
layer, but this remains a speculation. 

All the radio-occultation profiles displayed correspond to a uniform distribu- 
tion of the CH 4 mole fraction at 0.015. To check the sensitivity of the pro- 
files to variable CH 4 concentrations, temperatures were retrieved for the T14 
ingress (33° S) assuming the vertical variation obtained by the Huygens Gas 
Chromatograph Mass Spectrometer (GCMS): being « 0.015 down to 40 km, 
increasing to 0.024 at 20 km, then increasing more rapidly with decreasing 
altitude to 0.057 near the surface (Niemann et ah, 2010). The effect, illus- 
trated in panel / (blue dotted curve), is small. The differences are largest in 
the lowest 10 km, and the profile with variable CH 4 is about 0.6 K warmer 
at 5 km altitude. The retrieval with the GCMS CH 4 distribution in panel / is 
nearly coincident with the HASI profile. 

For a uniform CH 4 distribution, saturation will most easily occur in the 
tropopause region, because the temperatures there are lowest. For the temper- 
ature profiles near 30° S with 0.015 uniform CH 4 , the relative humidity over 
CH 4 ice is < 87%; at 53° S it is < 96%. Condensed CH 4 and N 2 are partially 
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miscible in the liquid and probably solid phases (Omar et al., 1962). If CH4 
and N2 are mixed in the ice phase, then the saturation vapor pressure of CH4 
would be less than over its pure solid phase and saturation might occur in the 
tropopause region. However, there are no available data on the properties of 
the vapor-solid phase equilibrium of CH 4 -N 2 mixtures, and what occurs on 
Titan is controversial. Tokano et al. (2006) have argued that the HASI tem- 
perature and GCMS CH4 gas profiles above the lowest few km are consistent 
with saturation over liquid CH 4 -N 2 mixtures, but above the freezing level near 
15 km, CH 4 vapor seems to be governed by the saturation pressure over CH 4 
ice. 


[Figure 10 about here.] 
[Figure 11 about here.] 
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Appendix A: Titan Tables 

This appendix presents two tables, Table 2 and Table 3. These tables contain 
numerical values for the temperature, pressure, and refractivity we obtain 
for the atmosphere of Titan for the occupations discussed in the main text 
at specific altitudes in the atmosphere. These numerical values at specific 
altitudes are obtained by linear interpolation from a smoothed version of our 
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computed results, which are not uniformly spaced in altitude. The smoothing 
was done to reduce interpolation artifacts on the coarse altitude grid presented 
in these tables. 


[Table 2 about here] 

[Table 3 about here] 

Appendix B: Nocturnal Inversion 

The three soundings near 30° S are close to dry adiabatic over much of the 
lowest 2 km. This is particularly evident in the T12 and T14 ingress soundings 
(Fig. 11). Moreover, both exhibit a stable inversion in the lowest ~200 m. The 
local solar time of both soundings is shortly before 5 A.M. This suggests that 
cooling of the surface drives the inversion. One can readily show that radiative 
cooling of the atmosphere just above the surface is too inefficient to establish 
a stable inversion on diurnal timescales. For a grey absorption the relaxation 
time r(m) to a sinusoidal perturbation of vertical wavenumber m is given by 
(Goody and Yung (1989), Ch. 10) 

1 Ank v dB 

r(m) pC p dT 

where p is atmospheric mass density, C p is the specific heat at constant pres- 
sure, k v is the volume absorption coefficient, with units of 1/length, T is 
temperature, and B = it~ 1 <JbT 4: is the Planck function integrated over all 
frequencies; Ob = 5.67 x 10 ~ 8 Wm~ 2 s~ x K ~ 4 is the Stefan-Boltzmann con- 
stant. Assuming radiative equilibrium and using the temperatures retrieved 
from the Voyager radio-occultation soundings, Samuelson (1983) estimated 
that k v /p ~ 7.6 x 10 ~ b m 2 kg~ l just above the surface. With 1/m ~ 200 m, 


n k v — i/ m w 

1 tan (— ) , 

m k v 



r(m) ~263 days, much longer than Titan’s 15.95-day diurnal period. 

However, a modest amount of turbulent mixing just above the surface, can 
account for the stable inversion over diurnal time scales. An eddy viscosity 
can be estimated as 

D = h 2 /t , 


where h is the vertical scale of the boundary layer, ~ 200 m, and t is the time 
scale of the diurnal cycle, 15.95 d/27r fa 2.2 x 10 5 s. This yields D ~ 0.1 m 2 s~ 1 , 
which is comparable to the global tropospheric value estimated by Flasar et al. 
(1981). For a 300-rn boundary layer, Tokano (Tokano et ah, 2006) estimated 
a value of D that was an order of magnitude smaller. However he related 
the eddy viscosity to the height of the Ekman layer where the surface winds 
are nearly geostrophic, which is a factor of 7 r larger than the vertical scale of 
the layer. Hence, for a given layer height and time scale, his estimate of D 
gave a value n' 2 smaller, which accounts for the difference. Note that the above 
estimate of D pertains to the stable inversion layer adjacent to the surface. The 
adiabatic layers above this likely have larger values, at least during daylight 
hours. 

Appendix C: The Effects of Ray Crossing from Sinusoidal Structure 
in an Isothermal Atmosphere 

The mean atmospheric structure typically has vertical profiles of density and 
refractivity that increase with decreasing altitude, on a scale comparable to the 
pressure scale height, H. Consider the simple case of an isothermal atmosphere 
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with constant H. When the bending angle, a , is small, 


a ~ 



No exp ( 



5 


where N 0 is the refractivity at z — 0, and R is the radius from the center of the 
occulting body to the ray periapsis point (Baum and Code, 1953). The domi- 
nant variation in the bending angle with altitude is from the exponential term, 
he., the bending angle to first order increases linearly with the refractivity. 

Now superpose a sinusoidally varying component to the exponentially vary- 
ing refractivity, with fractional amplitude e and vertical wavenumber m. The 
bending angle now does not necessarily increase monotonically with decreas- 
ing z. At altitudes where the perturbation increases with decreasing z, a will 
be larger than in the absence of the perturbation; where the perturbation de- 
creases with decreasing z, a will be smaller. This leads to the possibility of 
ray crossing. French and Lovelace (1983) found that crossing could occur in 
principle when 

(mH) 3/2 e > 1. 


Ray crossing becomes more likely for perturbations with larger amplitudes 
and smaller vertical scales. 

French and Lovelace (1983) were concerned with the ambiguities that can arise 
when inverting light curves from stellar occultations. When ray crossing oc- 
curred, their model light curves exhibited intense spikes. To illustrate the effect 
of ray crossing on retrieved temperature profiles when using an Abel inversion, 
we have used their algorithm and calculated the light curves in an exponential 
atmosphere for several cases with sinusoidal perturbations of different arnpli- 
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tudes. The light curves generated by ray tracing in the model atmospheres 
were inverted by means of an Abel transform. Figure 12 illustrates the input 
sinusoidal and retrieved temperatures, and also the vertical gradients of the 
retrieved temperatures. For low amplitudes, there is no ray crossing, and the 
Abel inversions faithfully reproduce the effects of the perturbed atmosphere 
structure. However, when the fractional amplitude becomes 0.03 or larger, dis- 
tortions become apparent. At e = 0.1, dT/dz exhibits a spiky structure that is 
reminiscent of that in Fig. 11. The vertical wavelength in Fig. 12 is H / 2 a; 10 
km, whereas the structure in Fig. 11 is smaller, for example, 300—400 m for the 
strongest spikes in the T12 ingress panel, about 30 times smaller. The above 
scaling suggests that the spiky structure observed in Titan’s lower atmosphere 
may correspond to fractional amplitudes of order 0.1/30 3 / 2 = 6 x 10~ 4 , or 0.06 
K. It should be noted that this corresponds to the average structure over a ray 
path distance comparable to the geometric mean of the vertical wavelength 
and Titan’s radius, several tens of kilometers. Nonetheless, small perturba- 
tions can have a noticeable effect on the Abel inversions, if the vertical scales 
are small enough. 


[Figure 12 about here. 
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Table 1: Titan Radio Occupations through 2006 



Ingress 

Egress 

Latitude 

Distance to 
Limb (km) 1 

Latitude 

Distance to 
Limb (km) 1 

Voyager 1 

12 November 1979 2 

6.2° N 

8700 

8.5° S 

19500 

Cassini 

T12: 19 March 2006 

31.4° S 

3745 

52.8° S 

5549 

T14: 20 May 2006 

32.7° S 

4027 

34.3° S 

4183 


1 Spacecraft distance to the limb corresponds to the time when the ray from the 
spacecraft to the receiving station on the Earth just grazes the solid surface of 
Titan, i.e., the “last” ray on ingress or the “first” ray on egress. Latitude is that 
corresponding to the periapsis of of the grazing ray. Dates correspond to Universal 
Time on reception. 

2 Voyager 1 results from Lindal et al. (1983) 


























Table 2: Values of the temperature (T, in K), pressure (P, in mbar) and refractivity (N) 
at selected altitude (h, in km) levels for the T12 ingress and egress soundings 








Table 3: Values of the temperature (T), pressure (P) and refractivity (N) at selected 

altitude (h) levels for the T14 ingress and egress occultations 
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Figures 



Figure 1. The power in the signal for the T12 ingress (a), T12 egress(b), T14 ingress 
(c), and T14 egress (d) soundings in the S, X, and Ka bands. The T12 power shown 
here was recorded at DSS-63 (S, X) and DSS-55 (Ka), and the T14 power was 
recorded at the same antennas. The black lines are a prediction of the power loss 
due to refractive defocusing computed as the refractivity profile was obtained. The 
S band power loss is essentially all due to refractive defocusing. Both X and Ka are 
visibly affected by N 2 absorption, with the Ka band power lost altogether within 
~ 10 km of the surface. The time axis is Earth received time in seconds past UT 
midnight on March 19, 2006 for T12 and seconds past UT midnight on May 20, 
2006 for T14. 






W. Longitude (°) W. Longitude (°) 

Figure 2. Latitude vs. longitude for each of the soundings discussed in this paper: 
(a) T12 ingress, (b) T12 egress (c) T14 ingress, (d) T14 egress. The (north) latitude 
and (west) longitude of the point on each ray closest to the surface of Titan (the 
ray periapsis) is shown as subsequent rays approach the surface of Titan (or in the 
case of egress soundings, recede from the surface). The circles denote the location 
of the indicated altitudes in km relative to a sphere of radius 2575 km. 
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Figure 3. Top. Temperature-altitude profiles for the T12 occultation, assuming a 
uniformly mixed atmosphere with mole fractions 0.985 N 2 and 0.015 CH 4 . The 
temperature was set to 185K at 300 km. Random errors (3<r) from thermal noise 
and spacecraft velocity errors are depicted on the right. The ephemeris was adjusted 
by using an along track shift so that the stratospheric profile agreed with profiles 
from CIRS (Achterberg et al., 2008b), in an attempt to minimize the effects of the 
velocity errors. The ingress ephemeris was shifted by -0.005 s, and the egress -0.014 
s. Both of these shifts produce velocity shifts outside the lcr quoted error for the 
ephemeris, although they are well within the 3 a level. Bottom. The same, except for 
T14. The shift for the ingress was -0.005 s, while that for the egress was -0.007. Both 
of these shift produce velocity shifts of ~ 1 mm/s, well within the quoted errors for 
the ephemeris. 



Figure 4. Top: The effect of starting the hydrostatic equilibrium computation at 
three different temperatures at 300 km. Note that the difference is small below 
~ 150 km but noticeable above that. Bottom: We start the integration at 148 km at 
172.5 K, approximately the same as the top of the HASI direct measurements. We 
compare our profile to the HASI direct measurements below 148 km. Both panels 
show the T14 egress occultation using 16 kHz X band data. 
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Figure 5. Doppler residual frequencies for DSS-63, X band, 16 kHz data, for the 
T14 ingress occultation. 




Figure 6. The ephemeris was shifted by the indicated times in seconds to produce 
these temperature profiles. This example is the T14 egress occultation using 16 kHz 
X band data. 



Figure 7. The results of a Monte Carlo study using the covariance matrix for the 
T14 egress occultation, using 16 kHz X band data. The velocity was shifted in a 
random direction by 3 times the covariance matrix elements. 100 cases were done. 
The solid line is the T14 egress profile, and the dotted lines encompass the 3 a error 
envelope. The maximum error is ±20K. 






Figure 8. Temperatures from the T12 and T14 occultations, obtained by adjusting 
the the spacecraft timing on ingress and egress until the retrieved temperature 
profiles in the upper stratosphere matched the CIRS zonal averages (see Fig. 3 
caption). Note that the vertical axis in this case is pressure, not altitude. 



Figure 9. The effects that random noise of magnitude similar to that already existing 
in the signal has on the temperature - altitude profiles, for T14 egress, 16 kHz X 
band data, DSS-63. The solid line is the T14 egress profile, and the dotted lines 
encompass the 3a error envelope. 








Figure 10. (a) Retrieved profiles of temperature from the T12 and T14 soundings for 
an assumed uniformly mixed atmosphere with mole fractions of N 2 : 0.986 and CH 4 : 
0.015. (b). T14 profiles compared to those retrieved from Voyager radio occultations 
(rescaled to an atmosphere with the same mole fraction of CH 4 ) and temperatures 
from the HASI experiment. The 1 -a errors indicated are from thermal noise only, 
(c), (d) Enlargement of (a), (b) in the troposphere and tropopause region, (e), (f) 
Enlargement of (c), (d) in the lowest 10 km of the troposphere. The retrieved T14 
ingress temperature assuming the variable profile of gaseous CH 4 obtained by the 
Huygens GCMS experiment is also shown in (f). The latitudes indicated for each 
occultation correspond to the periapsis of the rays grazing Titan’s surface. Dry 
adiabats are shown. 
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Figure 11. Vertical gradients of temperature in the lowest 10 km for the T12 and T14 
soundings. The vertical straight line in each panel corresponds to the dry adiabat. 
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Figure 12. Retrieved temperatures and their vertical gradients from Abel inver- 
sions for an isothermal atmosphere with a sinusoidal perturbation with fractional 
amplitude indicated in each panel. Temperatures and gradients are normalized by 
the mean isothermal temperature, Tq. Dotted curves. Input perturbation temper- 
atures, used in the forward ray-tracing calculation of light curves. Solid curves. 
Temperatures retrieved from the light curves using Abel inversions and the asso- 
ciated vertical gradients. Zero altitude corresponds to the half-power point on the 
light curve. The model calculations do not include a lower bounding surface. After 
French and Lovelace (1983). 









